Phase diagram and critical behavior of the three-state majority-vote model 
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The three-state majority-vote model with noise on Erdos-Renyi's random graphs has been studied. 
Using Monte Carlo simulations we obtain the phase diagram, along with the critical exponents. 
Exact results for limiting cases are presented, and shown to be in agreement with numerical values. 
We find that the critical noise q c is an increasing function of the mean connectivity z of the graph. 
The critical exponents f3/v, j/i> and 1/v are calculated for several values of connectivity. We also 
study the globally connected network, which corresponds to the mean-field limit z = N — 1 — ► 00. 
■ Our numerical results indicate that the correlation length scales with the number of nodes in the 

r ^ ' graph, which is consistent with an effective dimensionality equal to unity. 
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O ■ I. INTRODUCTION 



The majority- vote model with noise defined on a regular lattice is a system of spins where each one is allowed to 
be on two states only [H-0|- In this two-state model, each spin assumes the state of the majority of its neighboring 
spins with probability (1 — q) and the opposite state with probability q. The system presents an order-disorder 



phase transition as the noise parameter q reaches a critical value q c . Studies on the regular square lattice found 



q c = 0.075 ± 0.005 3j|, and critical exponents equal to those for the equilibrium Ising model in accordance with the 
conjecture by Grinstein et al. Q. 

The two-state majority- vote model (MV2) was also studied on a variety of complex networks [9j-H5j. such as 
undirected and directed random graphs 16], small- world networks [17j . and Barabasi- Albert scale- free networks [l8j . 
On undirected and directed random graphs, it was shown that q c is an increasing function of the connectivity of the 
graphs EH, [HI . On small- world networks, the critical noise is an increasing function of the rewiring probability 
[3, [l3[ . More general, these studies have shown that MV2 models defined on different complex networks belong to 
different universality classes and the calculated critical exponents depend on the topology of the complex network 
1^-21|. The generalization to a three-state majority vote model (MV3) on a regular square lattice was considered 
. [23 1 , where the authors found q c = 0.117 ± 0.001. The resulting critical exponents for this non-equilibrium 
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G\ ' MV3 model are in agreement with the ones for the equilibrium three-state Potts model [24]], again supporting the 
, conjecture of Ref. [8]. 

' In this paper we present an extensive study of the critical behavior of the three-state majority- vote model on Erdos- 
Renyi's random graphs [l6j . Monte Carlo (MC) simulations and standard finite-size scaling techniques are used to 
determine the critical noise parameter q c , as well as the exponents fi/v, "f/v and 1/v for several values of the mean 
connectivity z of the graph. We also study the globally connected network case. The phase diagram of the system is 
presented, and compared to our previously obtained diagram for the two-state model [l0| . Exact results for quantities 
of interest are obtained for the limiting cases q — > 0, and q — > 2/3, which agree with the simulation results. 

This work is organized in the following way: In section|TT]we describe the non-equilibrium three-state majority-vote 
model with noise, and introduce the relevant quantities used in our simulations. Sections IIIII and IIVI contain our 
jj] , results along with a discussion. Finally, in the last section we present our conclusions. 

a ■ 

II. MODEL AND FORMALISM 



The three-state majority- vote model with noise is defined by a set of spin variables {cj}, where each spin is associated 
to one vertex of an Erdos-Renyi's random graph and can have the values 1, 2, and 3. The connectivity of a vertex is 
defined as the total number of bonds connected to it, that is ki — J2j c iji where Cjj = 1 if there is a link between the 
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sites i and j and Cy = otherwise. A random graph is completely characterized by the mean number of connections 
per site, i.e. the average connectivity z, and the total number of sites N. 

The system evolves in time according to the following rules: For each spin we determine the state of the majority 
of its neighboring spins, that is, all the spins that are linked to it. With probability (1 — q) the new state of the 
spin agrees with the majority state of its neighbors and it disagrees with probability q, which is known as the noise 
parameter. In the case of a tie between the three possible states, each state is chosen with equal probability 1/3. 
In the case of a tie between two majority states, the spin assumes each one of these states with equal probability 
(1 — q)/2, and the minority state with probability q. Finally, in the case of a single majority state, the two minority 
states occur with equal probability q/2, and the majority state with probability (1 — q). It is clear that the rules just 
described present the C 3lJ symmetry with respect to the simultaneous change of all states a. 



Let k^ be the number of neighbors of site i in state a = 1, 2, 3, therefore k. 



(i) 



.(2) , lO) 



k\ — ki. According to the 



above rules we can write the following probabilities for a given spin to assume the state 1: 



P(l 
P(l 
P(l 
P(l 
P(l 



fcW = J fcC0 = fc(3) ) = 1/3 

k { p =kf ] >kf) = {l-q)/2 
k f) <k f) =} Jf) ) = q 

fcf) >fc f) )fc (3) ) = 1 _ g 

kP,k^ <kP)=q/2. 



(1) 



The probabilities for the other two states are obtained by the symmetry operations of the C% v group. For example, 
let us consider a neighborhood corresponding to the fourth and fifth rules, where we have a single majority state. In 
this case, the two minority states (say, states 2 and 3) occur with equal probability q/2, and the majority state with 
probability (1 - q). We can write P(l|fc(l) > fe(2),A;(3)) = 1 -q from the fourth rule, and P(2|fc(2), fc(3) < jfe(l)) = 
P(3|fc(3), k(2) < k(l)) = q/2, from the fifth rule. It is worth mentioned that the condition 1 — q (= Probability of 
choosing the majority state 1) > q/2 (= Probability of choosing a minority state, either 2 or 3) is valid for q < 2/3, and 
we conclude that q = 2/3 is the limit value for the noise parameter in the present three state MV model. Moreover, 
the probabilities defined by Eq.© satisfy P(l| . . . ) + P(2| . . . ) + P(3| . . . ) = 1. 

To study the critical behavior of the model we consider the magnetization Mjy, the susceptibility xn, and the 
Binder's fourth-order cumulant Un- These quantities are defined by 



M N (q) = «m) 



(2) 



XN(q)=N[{ (m 2 ) t ) c -((m) t ) 2 ] 



(3) 



U N (q) = 1 



( (m 4 ) t ) c 
3( (m?) t )l 



(4) 



where N is the number of vertices of the random graph with fixed z, (...) t denotes time averages taken in the stationary 
regime, and (...) c stands for configurational averages. In Eqs. ©-(j4]) m is defined in analogy to the magnetization 
in the three-state Potts model as the modulus of the magnetization vector, that is m = (m\ + to 2 + mj) 1 / 2 , whose 
components are given by 



(5) 



where the sum is over all sites in the graph, 8{a, ai) is the Kronecker delta function, and we introduce the factor 
\/3/2 in order to normalize the magnetization vector. 
In the critical region we assume the following finite-size scaling (FSS) relations [25[ 



M N (q) = N-P^MisN 1 ^) 



(6) 



XN (q) = N^xieN 1 ^) 



(7) 



U N (q) = UieN 1 /*) 



(8) 
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where s — q — q c , and the universal scaling functions M, \ and U only depend on the scaled variable x = eN 1 ^ . The 
above FSS relations follow from the ansatz that the correlation length scales with the number of nodes in the graph, 
that is £ ~ N, which is consistent with an effective dimensionality equal to unity. 

From the size dependence of Mm and xn we can obtain the exponents [3/9 and "f/9, respectively. The correlation 
length exponent 9 is calculated from the size dependence of the derivative of Binder's fourth-order cumulant with 
respect to the noise parameter, U' N (q — q c ). Furthermore, we use the hyperscaling relation 

2(3/9 + 1 /D = D eff , (9) 

to estimate the effective dimensionality of the system D e ff in order to check the FSS prediction for an effective 
dimensionality equal to unity. 

III. RESULTS 
A. Exact Results 

First we notice that there are only two independent components of the magnetization vector m a , since they obey 
the relation 

m 1 +m 2 +m 3 = 0. (10) 

Moreover, the norm of the magnetization is invariant with respect to any C 3v group symmetry operation. 

Let us consider the limit q — ¥ 0, where the probability of agreeing with the majority state equals unity. In this 
situation, after a transient period the system reaches the ordered steady state. Without loss of generality we can 
assume <Ji = l for all sites. In this case m\ — y/2/3 and m.2 = m 3 — — 1/\/6, thus Mjv(0) = 1. It is also straightforward 
to check that xjv(0) = 0, and U N (0) = 2/3. 

In the opposite limit (q — > 2/3) the system reaches a disordered steady state, where the average of each component 
of the magnetization vanishes. In fact, in this limit the probability of a given spin agreeing with the majority of its 
neighbors equals the probability of it agreeing with any of the other two minority states, i.e. (1 — q) — q/2. It is 
possible to write the probability distribution for the order parameter as a Gaussian distribution in the form [22| 

P(m) = -e- am \ (11) 

where a = l/(m 2 ). From this distribution we obtain that M N (q -> 2/3) = ^/n/2(m 2 ) 1/2 , XN(q -> 2/3) = N(4/ir - 
l)Mff, and U N (q -> 2/3) = 1/3. Moreover, it follows that (m 2 ) - N' 1 , M N (q -+ 2/3) - N^ 1 / 2 , and Xn(q -> 2/3) ~ 
N . These exact results as well as the predict size iV-dependence of the relevant quantities are in agreement with 
numerical results from simulations for all networks considered. 



B. Simulation 

We begin our simulations generating a random graph of size N and mean connectivity z in a disordered 
configuration where the state of each spin is 1,2 or 3 with the same probability. We used systems of size 
N = 1000,2000,4000,6000,10000,50000, and varied z from 1 to 50. To perform the dynamics we choose a site 
at random and, for a given fixed value of the noise parameter, we update its state in accordance with the dynamics 
rules given by Eq. (JTJ. A Monte Carlo step (MCS) is defined as N updates. We waited N r MCS needed for the 
system to reach the steady state, and the time averages, (...){, were estimated from the next N s MCS. The values 
of MCS used vary with N, z and q, typically we used N r > N s > 5000 MCS. For all sets of parameters, we have 
generated at least 100 distinct random networks in order to calculate the configurational averages (...)c- 

Fig. [T] shows the magnetization Mjv and the susceptibility Xn as functions of the noise parameter. The data were 
obtained from simulations on random graphs with N = 4000 sites and several values of the average connectivity z. 
In part (a) each curve for Mn, for a given value of z, clearly indicates that there exists a phase transition from an 
ordered state to a disordered state where the magnetization vanishes. We also notice that the transition occurs at 
a value of the critical noise parameter q c , which is an increasing function of the mean connectivity z of the random 
graph. In part (b) we show the corresponding behavior of the susceptibility xn- The value of q where xn has a 
maximum is here identified as q c (N). 

We also perform simulations for globally connected networks, that is, for random graphs with N nodes and con- 
nectivity z — N — 1. In the thermodynamical limit z = N — 1 — > oo, the magnetization is given bym = — |(q — |) 
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FIG. 1: Dependence of the magnetization (a) and the susceptibility (b) on the noise parameter q, for N = 4000 nodes. From 
left to right we have z = 2, 4, 6, 8, 10, 20, 30, and 50. In (a) the dashed line corresponds to the mean-field result m = — §(<?— §), 
in the thermodynamical limit z = N — 1 — > oo. In part (c) we plot the size dependence of the magnetization at q c — 2/3 for 
the case of z = N — 1, the globally connected network. 




FIG. 2: Binder's fourth-order cumulant as a function of q, for system sizes N = 1000, 2000, 4000. From left to right we have 
2 = 8, 10, 20, and 30. The horizontal line indicates the critical value U* = 0.42. The mean-field value is U* = 1/3. 

(the dashed line in Fig. 1(a)), from which we obtain the mean- field values qi MF ^ = | and /3 = 1. Fig. 1(c) shows 

the dependence of the mean-field magnetization at qi MF ' on the system size. The straight line confirms the scaling 
relation given by M^(q = 2/3) ~ iV -1 / 2 . The slope of the resulting straight line yields the exponent j3/v — 1/2. A 

similar analysis for the susceptibility at qi MF ^ yields Xn(q = 2/3) ~ N°, that is, j/v = 0. Note that the critical 
behavior for the globally connected network is in agreement with the exact results discussed above in the limit case 
of q = 2/3. The mean- field results for the critical noise parameter and critical exponents are given in Table |U 

In Fig. [2]we plot Binder's fourth-order cumulant Un for different system sizes ./V and four distinct values of z. The 
critical noise parameter q c , for a given value of z, is estimated as the point where the curves for different values of ./V 
intercept each other. We also obtained U* = 0.42(2) for the critical value of the cumulant at q c , which is independent 
of the connectivity z of the graph. The dependence of q c on z yields the phase diagram for the MV3 model shown in 

Fig. m 

The phase diagram of the MV3 model on random graphs shows that for a given graph (fixed z) the system becomes 
ordered for q < q c , whereas it has zero magnetization for q > q c . We notice that the increase of q c is more pronounced 
for small values of z. The error bars in q c (see Table HJ are much smaller than the symbols. In the figure, it is also 
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q c o.3 




FIG. 3: The phase diagram for the three-state majority-vote model (this work), showing the dependence of the critical noise 
parameter q c on the average connectivity z of the random graph. The phase diagram for the two-state model (Ref. [io( P is also 
included for comparison. 



included the corresponding phase diagram for the MV2 model obtained from Monte Carlo simulation in our previous 
work [IoT |. For both models the system exhibits an ordered state for all values of the mean connectivity greater than 
one. This is in agreement with the limiting value of z = 1 for the existence of a percolating cluster and, therefore, the 
onset of long-range order in the system. However, when z-}oowe obtain the upper limits q c = 0.5 and q c = 2/3, for 
the MV2 and MV3 models respectively. 

Fig. 0] shows the dependence of the derivative of the Binder's fourth-order cumulant at q = q c on system size. For 
clarity we have added in each curve the respective value of z. The straight lines, obtained from simulations with 
different values of the mean connectivity z, confirm the scaling relation given by Eq. ©. For fixed z, the slope of the 
resulting straight line equals the exponent 1/v. The results displayed in Table Q] indicate a weak dependence of the 
correlation length exponent with z. 

In order to obtain independent estimations for the critical noise parameter q c {z), as well as the exponent ratios j3 /v 
and j/v, for different values of the mean connectivity z, we consider the functions $ and \t defined as 

^NuN^-b-Hn^ (12) 



* Nl , N2 =b-Hn^ (13) 

where b = ln(iV2/iVi). The above functions relate the magnetizations and susceptibilities calculated with two different 
system sizes, N% and AT 2 . In fact, substituting the finite-size relation ([6]) into Eq. (IT21) . we obtain 

* NuN2 =P/D-b-Hn M ^ N \2 - (14) 
M(eNl /u ) 

At the critical value q c {z), the last term vanishes and we obtain Q(q c ) = ft I v. Analogously, substituting the finite-size 
relation ([7]) into Eq. (fT3|). we have ^{q c ) = l/v. 

In Fig. [5] we show the dependence with noise of the functions $ and "J, for z = 8 and several values of Ni and 
N2- From the intersection points in Fig. [3a) and Fig. [5jb) we obtain independent estimations for the critical noise 
parameter q c (z). For all networks considered in the simulations, we have obtained a quite satisfactory agreement 
between the two values of q c (z) determined in this way and the corresponding ones that follow from the analysis of 
Binder's cumulant (Fig. Moreover, we used the relations $(g c ) = @/v and ^(q c ) — j/v to calculate the exponent 
ratios for different values of the mean connectivity z. Table U shows the results for different values of the mean 
connectivity z. We call the readers attention to the difference between the calculated values of the exponents in the 
case of z — 2 and the corresponding results with z > 2. This might be an indication that, for z = 2, we need to take 
into account logarithmic corrections to the finite-size scaling relations [26| . 
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FIG. 4: (a) Plot of \n(U'(q c )) + z vs In AT. The exponent 1/v, for a given z, corresponds to the slope of the straight line obtained 
from a linear fit to the data. From bottom to top, z = 2,4, 6, 8, 10, 20, 30. 
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FIG. 5: The functions <3> and $ for the case of mean connectivity z = 8, and five different pairs of iVi and JV2- The intersection 
points give independent estimations for /3/v, j/v, and q c . The curves are cubic polynomial fitting to the data. 



Fig. [S] shows the data-collapse plot for M{x) = Mjq{q)N^/ v , which is a universal function of the combined variable 
x = N l ' v (q — q c ). We have also obtained quite good data-collapse for x(x) = \n {q)N~ 1 ^ v . The collapsing of curves 
for five different system sizes corroborates the quoted values for q c , /3/i>, j/u and 1/9. 

On Fig. |7| we present two different ways to obtain the data-collapsing for the universal scaling function U(x). Part 
(a) shows the standard data-collapse that follows from simulations with different values of system size N, for the case 
of mean connectivity z — 8 fixed. In part (b) we have fixed N = 4000 and used the data from simulations for varying 
connectivity. It is worth mentioning that in the last case the collapse for different values of z was obtained by using 
the values for the critical parameter q c and the exponent 1/v for the corresponding value of z (see Table UJ). In general, 
the universal scaling functions only depend on the scaled variable x — eN 1 ^ for a given system. However, the two 
data-collapsing in Fig. [3 indicate that U does not depend on the specific value of z, contrary to what is observed for 
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FIG. 6: The universal scaling function M versus the scaled variable x = eN 1 ^ . Data-collapsing for five different values of N, 
with z — 8. 
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FIG. 7: Data-collapsing of the Binder's cumulant. In (a) we show the collapse for five different values of N, with z 
In (b) the collapse is obtained for seven different values of z, with N = 4000 fixed. 
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IV. 



DISCUSSION 



Table U summarizes the values (along with errors) of the critical noise parameter q c , the critical exponents /3/9, 
j/p and l/P, and the effective dimensionality of the system. For all z considered, including the mean- field limit 
z = N — 1 — » oo, the value D e ff ~ 1 follows from the hyperscaling relation (Eq. (O), in agreement with the scaling 
ansatz for the magnetization and susceptibility Eqs. d6j [TJ) . In fact, since our original work [lfj several authors have 
studied different spin models on varied complex networks, always finding an effective dimensionality equal to one 
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TABLE I: The critical noise q c , the critical exponents, and the effective dimensionality D e ff, for the MV3 model on random 
networks with mean connectivity z. The mean-field (MF) results are shown in the last line. 



z 


Qc 


/3/P 


1 1 


l/v 

1 


D eff 


2 


0.084(5) 


0.15(3) 


0.75(2) 


0.81(6) 


1.05(8) 


4 


0.228(1) 


0.20(1) 


0.65(1) 


0.94(5) 


1.05(3) 


6 


0.3015(5) 


0.198(5) 


0.66(5) 


0.89(7) 


1.06(6) 


8 


0.3458(2) 


0.2105(5) 


0.62(5) 


0.895(5) 


1.04(9) 


10 


0.3785(2) 


0.22(2) 


0.68(2) 


0.90(2) 


1.12(6) 


20 


0.4586(5) 


0.205(5) 


0.67(2) 


0.95(1) 


1.07(3) 


30 


0.4957(1) 


0.22(1) 


0.66(2) 


0.95(2) 


1.1(3) 


50 


0.533(1) 


0.22(2) 


0.64(7) 


0.92(2) 


1.08(6) 


MF 


2/3 


1/2 





1/2 


1 



[ 1 ll — t 1 5I | . Even though this result seemed surprising at first, it is a direct consequence of the scaling for the correlation 
length £ ~ N 1 / D "ff , with D eff = 1. 

There are no previous works studying the three-state majority-vote model on Erdos-Renyi's graphs, to allow a 
direct comparison of the present results. Yet, for completeness, it would be of interest to mention earlier simulations 
of the majority-vote model on other kinds of networks. The only works on the MV3 model to this date, considered 
the model on a regular square lattice [22[ |23|. They find the universal value of Binder's fourth-order cumulant to be 
U* ~ 0.61, and the exponents (3/D = 0.134(5) and j/D = 1.74(2), all of which are in agreement with the results for the 
equilibrium three-state Potts model [24(. The present simulations of the MV3 model on globally connected networks 
yielded U* = 1/3, and the following mean-field exponents: f3/D — 1/2, 7/1/ = 0, /3 = 1, T> = 2, and 7 = 0. From our 
simulation results we can conclude that the MV3 model defined on a regular square lattice, on Erdos-Renyi's random 
graphs, and on the corresponding globally connected network (the mean-field limit) belong to different universality 
classes. 

Comparing the current results with the ones previously obtained for the two-state model on random graphs 
we first notice that the ordered region (q < q c ) in the phase diagram of Fig. [3] is larger for MV3 than for MV2. This 
is expected since now we might obtain three possible majority states. It should also be clear from the calculated 
exponents that these models do not belong to the same universality class. In particular, U* ~ 0.30 for MV2, and 
U* ~ 0.42 for MV3. 



V. CONCLUSION 



We have obtained the phase diagram and critical exponents of the three-state majority-vote model with noise on 
random graphs. The second-order phase transition which occurs in the model with mean connectivity z > 1 has 
exponents that show a slight variation along the critical line. Nevertheless, our Monte Carlo simulations provide an 
effective dimensionality -D e // equal to one for all values of z. This result, which is in agreement with several previous 
studies on spin models defined on complex networks, was shown to be a consequence of the ansatz that the correlation 
length scales with the number of nodes. Future work on the two- and three-state Potts model on random graphs 
would be of interest in order to provide a direct comparison with our results in light of the conjecture by Grinstein et 
aZ., which states that reversible and irreversible models with same symmetry belong to the same universality class. 
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